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Abstract. The aim of this theoretical and modelling paper is to derive knowledge on the global and structural parameters of 
low-mass stars using asteroseismology and taking advantage of the stellar collective behavior within open clusters. We build 
stellar models and compute the seismic signal expected from main sequence objects in the 0.8 - 1.6 Mq range. We first evaluate 
apparent magnitudes and oscillations-induced luminosity fluctuations expected in the Hyades, the Pleiades and the a Persei 
clusters. The closest cluster presents a feasible challenge to observational asteroseismology in the present and near future. The 
remainder of the work therefore focuses on the Hyades. 

We combine seismological and classical computations to address three questions: what can be inferred about 1) mass, 2) 
composition and 3) extension of outer convection zones of solar analogs in the Hyades. The first issue relies on the strong 
sensitivity of the large separation to mass. We show that seismic constraints provide masses to a precision level (0.05 Mq) that 
is competitive with the actual mass estimations from binary systems. Then large separations (Av) and second differences ( 62 V) 
are used to respectively constrain metal and helium fractions in the Hyades. When plotted for several masses, the relation of 
effective temperature (T cff ) vs large separation (Av) is found to be strongly dependent on the metal content. Besides this the 
second difference main modulation is related to the second ionization of helium. An accuracy in the helium mass fraction of 
0.02 to 0.01 can be achieved provided mass and age are accurately known, which is the case for a few Hyades binary systems. 
The second difference modulations are also partly due to the discontinuity in stellar stratification at the convective envelope / 
radiative core transition. They permit direct insight in the stellar structure. We compute acoustic radii of the convective bases 
for different values of the mixing length theory parameter q , mlt in convection modelling, i.e. different convective efficiency in 
the superadiabatic layers. For a given effective temperature we show that the acoustic radius changes with convection efficiency. 
This suggests that seismology can provide constraints on the extension of outer convection and also more generally on the direct 
approaches of convection and dynamical phenomena being currently developed. 

Key words. Stars : open clusters, asteroseismology. Stars : mass, composition, structure. 


1. Introduction 

Helioseismology achieved many successes in allowing di¬ 
rect determinations of the depth of the solar convection zone 
(Christensen-Dalsgaard, Gough & Thomson 1991) or photo- 
spheric mass fraction of helium (Basu & Antia 1995). In a 
comparable fashion the development of space-borne asteroseis- 
mic experiments within the coming years and the increasing 
number of ground-based observations should yield extremely 
precious data on other stars. In this paper we discuss some as¬ 
pects of the possible contribution of asteroseismology to the 
determination of mass, composition and structure of solar-like 
stars. We address the Hyades, Pleiades and a Persei open clus¬ 
ters. These clusters have been extensively studied from the¬ 
oretical and observational viewpoints (see Lebreton 2000 for 
a recent review) and their stars are among the targets of fu¬ 
ture asteroseismic missions such as EDDINGTON (Favata et 
al. 2000). Some of them will also be interesting targets for the 


recent and promising instruments HARPS (Pepe et al. 2003) 
or ESPaDOnS (Donati, Catala & Landstreet 2004). Contrary 
to field stars they offer the possibility to combine observations 
for stars having similar ages and compositions. In consequence 
they make it possible to exploit collective effects in seismic and 
stellar properties of stars with different masses. Finally they 
span a range of age and composition that should translate into 
significant differences in the seismic signals of solar-like stars. 

We compute stellar models and subsequently use them 
to predict oscillation frequencies. We use the CESAM code 
(Morel 1997) for pre-main sequence and main sequence mod¬ 
elling. Above 5800 K our version of the code relies on the 
OPAL2000 equation of state (Rogers, Swenson & Iglesias 
1996) and opacities (Iglesias & Rogers 1996). For lower tem¬ 
peratures the OPAL equation of state is replaced by the MHD 
equation of state (Mihalas, Dappen & Hummer 1988) and the 
opacities are from Alexander & Ferguson (1994). The atmo¬ 
sphere is restored using Hopf’s law (Mihalas 1978). It is con- 
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nected to the envelope at optical depth 10 where the diffusion 
approximation for radiative transfer becomes valid (Morel et al. 
1994). Microscopic diffusion is always taken into account fol¬ 
lowing the Michaud & Proffitt (1993) prescriptions. Finally the 
convection zones are assumed to be fully homogeneous, and 
modeled using the mixing-length theory (hereafter MLT) in a 
formalism close to that of Bohm-Vitense (1958). Unless men¬ 
tioned explicitly we consider q-mlt = 1-766 which is our Sun- 
calibrated value. No overshooting has been considered in con¬ 
vective cores or outer convective regions. Once the stellar struc¬ 
ture has been obtained we use the adiabatic pulsation package 
of Jprgen Christensen-Dalsgaard (1982) to compute oscillation 
frequencies 1 . As usual n and f. denote respectively the radial 
order and the degree of the mode whose frequency is written as 
Vl, n - 

Before we present the approaches we used to relate the stel¬ 
lar parameters to the seismic signal we address the question of 
the modes amplitudes in the following section. This question is 
a crucial one for the efficiency of future seismic campaigns. 

In the Sect. 3 we focus on the impact of global parame¬ 
ters on the asteroseismic data. Both initial mass and chemical 
composition are important in stellar evolution. Helium is the 
second most abundant constituent of the Universe. It is is how¬ 
ever impossible to determine the helium fraction from spectral 
analysis of low-mass stars. Its abundance therefore remains un¬ 
determined in all open clusters except the Hyades. In this open 
cluster lower-main sequence fits that constrain helium content 
have been performed (Perryman et al. 1998; Pinsonneault et 
al. 1998). Stellar mass determination based on close binaries is 
similarly a difficult task. Except in the case of nearby systems 
such as a Cen the uncertainty in the mass generally is /, 0.1 Mq 
(Soderhjelm 1999). Addressing the specific case of the Hyades 
cluster we show what improvements seismology could give re¬ 
garding the helium mass fraction and more generally composi¬ 
tion. 

In Sect. 4 we discuss the seismic characterisation of the 
transition from the radiative interior to the convective envelope 
in solar-type stars. Then we evaluate the impact of changes in 
the convective efficiency in the superadiabatic regime over the 
depth of the convection zone. We therefore propose seismol¬ 
ogy as a tool to calibrate the ctmlt parameter or more impor¬ 
tantly the convection efficiency to be inferred from direct hy¬ 
drodynamic computations for different surface conditions (r e ff, 
log g). 

A discussion of our results and our conclusions are given in 
Sect. 5. 

2. Detectability of oscillations 

The detection and study of solar-like oscillations modes are re¬ 
stricted by several factors. Not only are they limited to a certain 
range of visual apparent magnitude but they also depend on 
the amplitude of the luminosity modulations. In this section we 
discuss the anticipated amplitude of the p-modes for the vari¬ 
ous stars we modeled. This is a central point in the perspective 

1 Aarhus Adiabatic Pulsation Package, presently available at 

http://astro.phys.au.dk/~jcd/adipack.n/ 


of space asteroseismic missions such as COROT (Baglin et al. 
2002) and MOST (Walker et al. 2003) or ground based obser¬ 
vations. We have derived the apparent magnitudes in V band on 
the ~ 0.77 to 1.6 Mq stellar mass range and for the three open 
clusters (TableQ. We take M Bo i = 4.72 - 2.5 log 10 (L/Lq) (ie 
M Bo1 q = 4.72) and we apply the bolometric correction BCv 
as deduced from a bilinear interpolation in T e g and log g tables 
provided in Houdashelt et al. (2000) 2 where T e ff is the effective 
temperature and log g the log of the surface gravity. We have 
used distance moduli based on HIPPARCOS parallaxes have 
been used. We take m—M = 5.37 + 0.07 and 6.31+0.08 respec¬ 
tively for the Pleiades and the a Persei cluster (van Leeuwen 
1999). For the Hyades we assume m - M = 3.33 ± 0.01 
(Perryman et al. 1998). The extinction A v is negligible for the 
Hyades, close to 0.04 mag for the Pleiades and 0.1 mag for a 
Persei (Pinsonneault et al. 1998). Owing to these values no red¬ 
dening has been taken into account in the process of computing 
visual apparent magnitudes from the absolute bolometric mag¬ 
nitudes provided by the models. The magnitudes we report are 
therefore only accurate to 0.1 mag. 

Let us now focus on the expected amplitudes for the modes. 
These are determined by the relative efficiencies of the driv¬ 
ing and damping processes. For solar analogs the most prob¬ 
able explanation of driving of modes consist of the fluctua¬ 
tions of the Reynolds stress associated with convection in the 
superadiabatic layers. In stellar conditions the acoustic emis¬ 
sion increases as M, 7 - 5 where M, is the turbulent Mach num¬ 
ber (Goldreich & Kumar 1990). Recent computations suggest 
that M, reaches a maximum around 1.6Mq on the zero age 
main sequence (ZAMS) and so does the oscillation amplitude 
(Houdek et al. 1999). The most massive stars we address here 
are therefore expected to have the largest fluctuation ampli¬ 
tudes. In terms of surface velocity, the observations suggest: 

ttosc OC (L*/M*) 0 ' 7 (1) 

In this relation, as in the following ones, L* and M* stand 
for stellar bolometric luminosity and mass in solar units respec¬ 
tively. Kjeldsen & Bedding (1995) and Houdek et al. (1999) 
have proposed similar laws with somewhat higher exponents. 
However they predict large velocity fluctuations that exceed the 
actual observations in the case of most massive stars such as 
Procyon (1.46 Mq). Using theoretical models of stochastic ex¬ 
citation (Samadi & Goupil 2001), Samadi et al. (2004) confirm 
this power-law dependence with an almost identical exponent 
(0.8 instead of 0.7). Space-borne experiments will measure 
brightness modulations into which v osc has to be converted. 
Such an analysis is performed in Kjeldsen & Bedding (1995) 
and we will not reproduce it here. As the main result of the 
computations one obtains {j^) B , K which together with 
Eq.Q yields an estimate of the dependence of the amplitude of 
the luminosity fluctuations on fundamental stellar parameters : 

(—) oc—] (2) 

l L JBoi V7W [M+l 

2 IAU Commission 25 (1997) XIII General Assembly recom¬ 
mended the bolometric magnitude of the Sun to be defined as +4.75. 
However, following Houdashelt et al. (2000) we use +4.72 for consis¬ 
tency. 
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Table 1. Apparent visual magnitude my and expected luminosity fluctuations for 0.78, 0.8, 0.95, 1, 1.15, 1.3, 1.5 and 1.6 Mq 
stellar models. Assumed ages are 90, 125 and 625 Myr for a Persei, Pleiades and Hyades respectively. We consider non-solar 
repartition of metals in the Hyades and Y=0.26 (see Sect.0. For the two other clusters we assume solar composition. Luminosity 
fluctuations are probably underestimated for the 1.5 and 1.6 Mq models. Rows 4 to 7 list various seismological features for the 
age and composition of the Hyades. 


Mass Mq 

0.78 

0.8 

0.95 

1 

1.15 

1.3 

1.5 

1.6 

my (Pleiades) 

12.1 

12.0 

10.9 

10.6 

9.8 

9.1 

8.4 

8.1 

m v ( a Persei) 

13.1 

12.9 

11.8 

11.5 

10.7 

10.1 

9.3 

9.1 

my (Hyades) 

10.1 

9.9 

8.8 

8.5 

7.7 

7.0 

6.3 

6.0 

Luminosity fluctuations (ppm) 

1.8 

1.9 

2.9 

3.2 

4.6 

6.2 

8.7 

9.9 

Frequency of maximum amplitude y max (mHz) 

3.5 

3.5 

3.9 

3.9 

3.7 

3.3 

2.7 

0.8 

Cut-off frequency v c (mHz) 

8.7 

8.5 

6.9 

6.3 

4.8 

3.6 

2.4 

2.1 


If we now use the Sun as a calibrator, this simple law pro¬ 
vides an estimation of the luminosity variations induced by 
oscillations. Such an assumption should be adequate, for the 
stars we consider have solar-like (convectively powered) os¬ 
cillations and the external convection properties should vary 
smoothly with stellar type. The oscillations associated with so¬ 
lar luminosity fluctuations are of the order of 4 ppm (Kjeldsen 
& Bedding 1995). Therefore if we express the stellar mass and 
luminosity in solar units we finally deduce : 



Using this last equation, Table^gives the amplitude of the 
luminosity fluctuations as a function of mass for the Hyades 
models. It also indicates the magnitudes of the corresponding 
objects in the other two clusters. Because of differences in com¬ 
position and evolutionary stage, the relative amplitudes of the 
fluctuations decrease slightly at a given mass when going from 
Hyades to the younger clusters. We do not compute any metal- 
licity impact on the amplitudes. The small variations are in¬ 
stead most probably due to an age effect which is expected to 
be very small for ~ 0.5 Gyr difference (Houdek et al. 1999) and 
we do not reproduce them here in detail. The maximum varia¬ 
tion concerns the most massive star where it reaches 0.6 ppm 
i.e. w 6 %. 

Also of interest is the range of frequencies for which solar¬ 
like p-modes oscillations should be the best visible. This fre¬ 
quency range roughly covers a region from the frequency v max 
where the largest amplitude is expected up to the cutoff fre¬ 
quency y c and a similar frequency extent below v max . We pro¬ 
vide v max in Tabled together with the atmospheric cut-off fre¬ 
quency v c ; v c is computed from its usual expression 27rv c = 
with c s and H p respectively the sound speed and pressure 
scale height in the atmosphere ; v max is estimated following 
Vmax = ^rr~ where Li max stands for the maximum convective ve- 

n p 

locity in the superadiabatic layers and H p is the corresponding 
pressure scale height. Because of the sharp drop in density and 
temperature on the one hand and the requirement for nearly 
constant convective flux on the other hand, convective motions 
are substantially accelerated in the upper superadiabatic layers. 


The driving of solar-like oscillations is highly sensitive to the 
kinetic energies of convective eddies. Thus its efficacy reaches 
a maximum in the upper superadiabatic layers. Modes of fre¬ 
quencies y max similar to the local turn-over frequency y to where 
the convective velocity is maximal exhibit the largest oscilla¬ 
tion amplitudes and therefore y max ~ y t0 . For modes having 
v 0 sc > v to the excitation vanishes because of a lack of coher¬ 
ence between eddies at higher frequencies; while in contrast for 
v 0 sc < v to there is a rough equipartition between convective ki¬ 
netic and mode energies (Goldreich & Keeley 1977). However 
because of the increasing of the inertia of the modes the am¬ 
plitudes diminish with the decrease of v osc . We made a rough 
estimate for the convective velocity v using the mixing length 
theory formulae (precise demonstrations are to be found in text 
books): 


v = 

B = 


O'MLT 


g$H p 


r 

W 2 


1 (a M LTH p ) 4 (pC p ) 2 6g 


162 


K 2 H n 


(4) 


where the various symbols have their usual meanings : ctmlt 
is the mixing length parameter; g is the local gravity; 6 = 
- F is the convective efficiency; C p is the gas spe¬ 

cific heat at constant pressure ; p and K are respectively the 
gas density and thermal conductivity. These various quantities 
are derived from our models. We use the Eq.0as a scaling law 
considering that in the case of the Sun the maximum oscillation 
amplitude is at 3.3 mHz which corresponds to the well known 
5-minute oscillations. Our solar reference model for convec¬ 
tion velocities is the calibrated ’Btz’ model from Brun, Turck- 
Chieze & Zahn (1999) for which the equation 0] provides a 
maximal velocity v max = 2.6 km s~ 1 . 

The magnitudes and oscillation amplitudes computed here 
are indicative of the observational accuracy needed to address 
solar objects in nearby open clusters using seismology. Even 
though the COROT space mission does not primarily aim at 
clusters we will consider its performance as indicative of what 
is now within technical reach. For a my =9 magnitude the 
threshold of oscillation detection expected for this mission is 
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{^) B j = 2.4 ppm (Baglin et al. 2001). Table 0 therefore sug¬ 
gests that Hyades stars around and above solar mass should 
become accessible to seismic measurements in the near future. 
A Hyades star of 0.95Mq is for instance interesting in several 
respects. 

It should exhibit my in the 8.7 to 9.1 range depending on 
the actual age and composition of the Hyades (see Table 0 of 
Section0. This somewhat high magnitude is however compen¬ 
sated by a (x) B i we ex P ect t° be above 2.8 ppm. In the next 
section we develop the informations on mass one could extract 
from such a seismic target. For the Pleiades and a Persei the 
situation is far less favourable because of their distances. With 
1 3Mq and 1.6 Mq stars having a visual magnitude of ~ 9 in 
these respective clusters, solar mass stars will be difficult tar¬ 
gets there. The oscillation amplitudes expected above 1.3Mq 
are sensibly higher than solar ones and render corresponding 
stars better targets. One nevertheless needs to be cautious con¬ 
cerning the granulation noise which may limit the observations 
of luminosity fluctuations. We end this section with two re¬ 
marks: 1) Because of its youth a Hyades solar analog has a 
slightly larger magnetic and chromospheric activity than our 
Sun (see Paulson et al. 2002 for recent evaluations). This could 
be problematic for the clear detection and and interpretation of 
the oscillatory modes. 2) Stars having masses between 1.5 and 
1.6 Mq are 5 Scuti variables when they are near to the ZAMS. 
We are presently not interested in the difficult question of low 
frequency modes excited by radiative processes. In order to 
avoid such modes we shall limit ourselves in subsequent com¬ 
putations to modes of orders above n=10 for which the k mech¬ 
anism has been shown to be inefficient (Michel et al. 1999). We 
note that both radiatively and convectively driven oscillations 
are expected in 6 Scuti variables. Furthermore because the re¬ 
gions of 6 Scuti involved in the driving of the modes are not at 
dynamical equilibrium the expected luminosity fluctuations ^ 
are much higher than predicted when using equation 0 (see 
Samadi et al. 2002). 

3. Global parameters and seismology in the 
Hyades 

3.1. Composition determination 

Asteroseismology may complement spectroscopic chemical 
determinations. It will provide some indications on composi¬ 
tion as, because of opacity effects, metals and helium have a 
direct impact on stellar global parameters (luminosity, effec¬ 
tive temperature) and structure. This information is of course 
dependent on the various assumptions made in the modelling 
(convection parameter, diffusion processes, age...) 3 . 

As a first step, we will exploit the collective seismic be¬ 
haviours of stars of similar age and composition. In open clus¬ 
ters it is indeed possible to draw conclusions from seismic HR 
diagrams in the same way as from the classical HR diagram. 
Figure 0 is a transposition of the usual theoretical HR dia¬ 
gram where bolometric luminosity has been replaced by the 

3 The past year’s developments concerning solar surface oxygen 
fraction however reminds us that spectroscopic measurements also are 
strongly dependent on models (Asplund et al. 2004). 


Table 2. The five composition considered for the Hyades. 



Cl 

C 2 

C3 

C4 

C5 

Helium mass fraction 

0.26 

0.25 

0.24 

0.28 

0.26 

[Fe/Hf (dex) 

0.127 

0.127 

0.127 

0.127 

0.127 

[O/Hf (dex) 

-0.07 

-0.07 

-0.07 

-0.07 

0.127 


a [Fe/H] represents iron and iron-group metals while [O/H] repre¬ 
sents oxygen and non-iron-group metals. 



6500 6000 5500 5000 

Effective temperature (K) 

Fig. 1. Large separation as function of effective temperature be¬ 
tween ~ 5000 and 6500 K for varying possible composition 
model of the Hyades. The solid line and dashed line corre¬ 
spond to Cl ([Fe/H] ± [O/H], Y = 0.26) and C5 ([Fe/H] = 
[O/H], Y = 0.26) respectively. For the same value of the ab¬ 
scissa they are separated by ~ 15 /./Hz. The dot-dashed line and 
the triple dot-dashed line represent C2 ([Fe/H] + [O/H], Y = 
0.25) and C3 ([Fe/H] =£ [O/H], Y = 0.24) respectively. These 
models show that similarly to the usual HR diagram a decrease 
of the helium fraction has the same effect as the increase of 
the metal fraction. The diamonds, crosses and triangles are the 
positions of the 0.95, 1.07 and 1.25 Mq models for each com¬ 
position. The horizontal error bar below the curves at ~ 5600 K 
shows the typical uncertainty of the effective temperature. We 
also indicated a 1 //Hz wide vertical bar. The expected preci¬ 
sion for mode frequency determination from EDDINGTON is 
0.3 //Hz. 

large separation. This figure therefore does not show the tradi¬ 
tional seismic HR diagram in which the small separation (de¬ 
fined as 6v{ y „ = ve.n - V{+ 2 ,n-\) is plotted as a function of the 
large separation (defined as Avf „ = Vf_ n - Vf M . \) (Christensen- 
Dalsgaard 1988). In this study we have considered five plausi¬ 
ble compositions for the Hyades, summarized in Table 0 The 
first composition (hereafter C1) assumes that [Fe/H] = 0.127 
dex (Boesgaard & Friel 1990) applies to all iron-group met¬ 
als we consider 4 but uses [O/H] = -0.07 (Garcia Lopez et 
al. 1993) for oxygen and other metals. It therefore presents a 

4 Ni, Fe, Mn, Cr, Ti 
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non-solar-fraction repartition among metals. For this Cl com¬ 
position we choose a helium fraction Y = 0.26 as supported by 
fitting of the lower main sequence (Perryman et al. 1998) and 
by mass-luminosity relation of binaries (Lebreton et al. 2001). 
The second (C2), third (C3), and fourth composition (C4) only 
differ from Cl through Y=0.25 Y=0.24 and Y=0.28 respec¬ 
tively. The fifth composition (C5) assumes [Fe/H] = 0.127 
dex, Y=0.26 and (the usual) solar repartition among metals. 
Provided that the effective temperature is reliably determined, 
Fig.^illustrates that it should be possible to make a clear dis¬ 
tinction between Cl and C5 : in the 5000 to 6000 K effective 
temperature range the mean large separation Av (see Sect. E3 
for the precise definition) varies by ~ 10//Hz between a metal- 
rich Hyades star model (considering that the metal fractions 
follow the [Fe/H] =0.13 dex value) and its metal-poor coun¬ 
terpart (considering [O/H]=-0.07 dex as Garcia Lopez et al. 
(1993) suggests) that exhibits the same effective temperature. 
Such a difference is therefore quite significant when translated 
to expected accuracies of seismic experiments. The typical 
uncertainty in effective temperature found in the Hyades low- 
mass stars is 55 K (Thorburn et al. 1993) whereas the large 
separation determination should be affected by an error around 
0.1 /.(Hz. The temperature uncertainties therefore have a much 
larger impact on the interpretation of Fig. [2 in terms of compo¬ 
sition of the Hyades. 


It is worth estimating the limits of this method. Although 
the main sequences in Fig. [3 strongly depend on composition 
caution is needed when interpreting their positions in these 
terms. The large separation is directly fixed by the radius and 
consequently by the assumptions about convection modeling. 
For instance for a Hyades Cl composition and 0.95 Mq model 
a change of umlt from our standard value 1.76 to 1.56 induces 
decreases of 70 K and 4 //Hz in T e ff and mean Av respectively. 
These variations of convective efficacy remain plausible (see 
Sect.^J. Let us now compare the accuracy of the determination 
of the composition from our seismic HR diagram with the pre¬ 
cision of the classical HR diagram. The smaller the uncertain¬ 
ties in the large separation (respectively the luminosity) the bet¬ 
ter defined the main sequence and the more precise the compo¬ 
sition determination through a fit of the lower main sequence or 
the age determination from the turn-off position. In the frame¬ 
work of usual HR diagram analysis, Stello & Nissen (2001) es¬ 
tablish that the mean error in absolute visual magnitude can be 

7 £Ylf 1/2 


expressed as cr(My) = 


no-(b-y) 1 + (2.\1- 


, where b 


and y are the color bands from the Stromgren ubvy photome¬ 
try while n is the parallax. The precision in My for individual 
members of clusters is mainly affected by uncertainties in these 
quantities. In the case of the Pleiades and the a Persei clusters 
the HIPPARCOS value for ^ of the brightest stars exceeds 
0.1 (Eggen 1998). On the other hand Stello & Nissen (2001) 
estimate cr(b-y) ~ 0.01. The global error in the absolute mag¬ 
nitude My of individual stars in the Pleiades and a Persei there¬ 
fore is i (0.12 2 + 0.21 2 ) 1/2 0.25. In Fig.|5]we illustrate and 

compare the impact of the composition on positions of the main 
sequence in a classical and in the seismic HR diagram. We con¬ 
sider three plausible compositions for the Pleiades that present 
variations of AT = 0.1 and A [Fe/H] = 0.05 dex among each 


other. We remark that with respect to the errors the separation 
along the ordinate axis is much larger in the seismic than in the 
classical HR diagram. 

Recently de Bruijne et al. (2001) derived extremely accu¬ 
rate individual Hyades distances making use of HIPPARCOS 
secular parallaxes. This brings the accuracy in the absolute 
magnitudes My down to the ~ 0.05 level and allowed the au¬ 
thors to produce an unprecedentedly narrow main sequence. 
My is however not a quantity provided through stellar structure 
computation and that can be used in a theoretical HR diagram. 
We note that the B - V to 7’ c fr calibration is still controversial 
(Bessell et al. 1998, Lejeune et al. 1998) and is found to be in¬ 
appropriate below ~ 5000 K (de Bruijne et al. 2001). Different 
transformations lead to zero age main sequences varying by up 
to 0.2 mag in the My, B-V plane and between B-V = 0.7-0.8 
(ie r eff ~ [5500,5200] K) whereas they are similar in the 
MboI, T e g plane (Robichon et al. 1999). Moreover the My trans¬ 
formation to bolometric luminosity is subject to the difficult 
question of bolometric corrections ( BCy). BCy depends on the 
effective temperature and also critically depends on metallicity 
(Alonso et al. 1996) and surface gravity (Arribas & Martinez- 
Roger 1988). We consequently stress here that as long as de¬ 
tailed metal abundances remain uncertain absolute luminosity 
determination might be biased. In view of this different sources 
of uncertainty we consider that seismic diagrams could give 
valuable insight into composition even for the nearest open 
cluster. 

In Figs. 0 and |3 the main sequences are not surprisingly 
dependent on the helium fraction. Slight opposite variations 
in the helium or metal fractions have similar impact on stel¬ 
lar luminosity and effective temperature so that a variation of 
the helium fraction or metal fraction do not seem to be a pri¬ 
ori distinguishable. Hence main sequence fitting can only lead 
to metal abundances once the helium fraction is known (and 
conversely). We nevertheless point out that if the Hyades metal 
fraction is around -0.07 dex instead of 0.1 dex then the main 
sequence clearly lies above this last metal-rich main sequence 
for any possible helium fraction : the C3 composition main se¬ 
quence is in a sense a limit that can be assigned to metal-poor 
main sequences for it has Y=0.24 which is critically close to the 
actual primordial helium evaluations : Izotov & Thuan (1998) 
propose T p = 0.245 + 0.004 and Peimbert & Peimbert (2000) 
propose T p = 0.2345 ± 0.003. 

In order to specifically check for the helium effects we have 
used the second difference : (Lvo = ve, n -i - 2v/> + v/>+i. As 
is well known, the fluctuations of 52Ve,n with frequency out¬ 
line the presence of two rapid variations in the stellar struc¬ 
ture : firstly the region of the second ionization of helium, sec¬ 
ondly the discontinuity at the base of the outer convection zone. 
In a recent work Basu et al. (2004) suggest that it is possi¬ 
ble to determine the helium fraction Y using the mean ampli¬ 
tude of the second difference modulation provided that both 
mass and radius are accurately known. In that work the signa¬ 
ture appears more clearly around 1.2 Mq (see also Lopes et al. 
1997). Masses are generally not known in the Hyades and there 
are no measurements of radii. However the age is known and 
some binary systems such as Finsen 342 have their masses pre¬ 
cisely determined. Because of its mass (1.25Mq) and luminos- 
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Fig. 2. Left panel: usual theoretical HR diagram for 
Pleiades. Continuous line: solar composition stars (X=0.7354, 
Y=0.2465); Dashed line: helium fraction is increased by 0.1 
but [Z/H] is kept unchanged so that X=0.7256, Y=0.2565. 
Dotted line: [Z/H] =0.05 dex with respect to solar composition 
so that X=0.7275, Y=0.2524. The vertical error bar (0.2 mag) 
is representative of the absolute magnitude uncertainty for the 
Pleiades cluster. Triangles, diamonds, and crosses show 0.8, 
0.9 and 1 M© stars respectively. Right panel: seismic HR di¬ 
agram for the same cluster: Large separation vs effective tem¬ 
perature. The plot conventions are similar. The minute vertical 
error bar of 1 //Hz is representative of the maximal uncertainty 
in the large separation. We note that the separation between 
main sequences is much larger in the seismic diagram. 


ity (my ~ 7.4) Finsen 342 B should be an interesting target for 
helium determination. We have built models of this star 5 and 
computed the second differences for different initial helium Y 
values. Assuming thar the second ionization of helium acts as 
a discontinuity, the second difference was fitted to the form 6 : 


62 V = + a 2 v + 


+ 



sin(47rvr Hf , + <p He ) 


+ ^ci + sin(47rvT BCZ + 4>bcz) (5) 

5 Interestingly the 9 1 Tau B star has a very similar mass (1.28 ± 
0.13 Mq) to that of Finsen 342 B so that our analysis also applies to 
this object. 

6 The region of the second ionization of helium has a non-negligible 
extent and one should use other fits than Eq.[5] However these fits 
provide similar helium fraction estimations as the one based on Eq.[5] 
see Basu et al. (2004) 
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Fig. 3. The mean amplitudes of the second difference modula¬ 
tion induced through the second ionization of helium scaled by 
Av/Avq. The means are given as function of helium fraction in 
the envelope at the age of the Hyades : Y=0.213, 0.229, 0.236, 
0.242, 0.248, 0.235 and for composition Cl (crosses). Because 
of the diffusion process these values do not correspond to the 
initial helium fractions that are indicated next to the crosses. 
Error bars of 0.1 //Hz are only indicative. 


Finally the amplitude A Ht , of the helium-induced modulation 
was taken as the mean value of b\ + over the the range 
[1.5,3]mHz, both quantities being scaled by Av/Av© (we take 
Av© = 135/vHz). This last factor is introduced because the 
amplitude of the signal scales with the mean density, as does 
the large separation Av (see Basu et al. 2004). Figure 0illus- 
trates the results: a variation in helium fraction from 0.24 to 
0.26 induces a change of 0.23/vHz in the fitted mean amplitude 
A H( ,. Let us make here an important remark. In Hyades stars 
of w 1.2M©, diffusion is definitely not a negligible process be¬ 
cause of the shallower outer convection zone. In the models at 
625 Myr, the initial outer convective region helium fractions 
of 0.24 and 0.26 have decreased to 0.229 and 0.242 respec¬ 
tively. The A Ane variation is therefore larger than 0. 1 /jHz for 
AT = 0.01. The diffusion effect becomes even more dramatic 
above T = 0.28 as the extension of the convection zone contin¬ 
uously diminishes with increasing helium fraction. The model 
that has an initial helium fraction of 0.3 indeed exhibits a lower 
Hyades age content (T = 0.235) in its convection region as the 
model having initial helium fraction of 0.28. In this instance we 
recall that the evolution of the A Ht , amplitude is obviously not 
only related to the helium fraction but also depends on physical 
conditions in the region of helium second ionization. 

3.2. Mass determination 

Binary systems have been extensively used to compute stel¬ 
lar masses on the basis of astrometric and spectroscopic mea¬ 
surements. A recent work of Soderhjehm (1999) based on 
Hipparcos astrometric data evaluates individual masses within 
27 systems. For the subsample of 14 systems within the ap¬ 
proximate mass range from 0.5 to 1.5 M© the mean precision 
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Fig. 4. Second difference and fit for Finsen342 B if a helium 
fraction of Y=0.28 is assumed. 

in mass determination is 0.11 Mq and exceeds 0.06 for all 
of the stars but two. In the Hyades cluster no system has its 
mass determined to better than 0.16 Mq and individual masses 
are available for only one system with a precision better than 
0.2 Mq. Making simultaneous use of spectroscopy and speckle 
interferometry can significantly improve the accuracy of such 
mass determinations : Pourbaix (2000) determines the masses 
to better than 5 percent for 25 stars from 40 binary systems. In 
very favourable cases the mass uncertainty has been decreased 
to slightly below O.OIMq (Latham et al. 1996). Combining ob¬ 
servations of the Hyades binary system Finsen 342 (70 Tauri) 
over a period of 15 yr, Torres, Stefanik & Latham (1997B) 
found 1.363 + 0.073 Mq and 1.253 + 0.075 Mq for the pri¬ 
mary and the secondary components. However, because of 
the observational difficulties and the long amount of time re¬ 
quired for careful interferometric and spectroscopic measure¬ 
ments only very few other stars in the Hyades (vB22, 51 Tau, 
6 2 Tau) have had their masses determined with similar preci¬ 
sion (Torres Stefanik & Latham (1997A); Peterson & Solensky 
1988). 

Asteroseismology opens up new possibilities for accurate 
absolute or relative mass determinations. Its methods are not 
restricted to close binary systems but extend the determina¬ 
tion of masses to all stars. We note by the way that the non¬ 
occurrence of oscillations could be used as an indication of 
possibly unresolved multiple systems that are sometimes sus¬ 
pected (Soderhjelm 1999) in classical observations. In the fol¬ 
lowing we focus on Hyades stars in the 0.77 to 1 Mq mass 
range. 

The asymptotic theory of p-mode frequencies (Tassoul 
1980) makes it possible to relate the large separation Aiq„ = 
ve,n - ve, n - 1 to the stellar acoustic radius r«, namely : 

1 r R dr 

T7- — = t r (6) 

2Av 0 , Jo c 

Table0shows the mean Av for masses varying from ~ 0.77 
to IMq and for possible ages and compositions of the Hyades. 
We compute Av as the average of Av/ : n over t = 0, 1 and 2 and 


Table 3. Stellar parameters for 0.769(vB22B), 0.8, 0.95 and 
1 Mq Hyades stars of composition Cl. For 0.95 Mq results 
are given for various ages and compositions. Large separation 
computations are based on the mean of Avc, n in the frequency 
range going 3500 to 4500 pHz. The scatter in this mean always 
lies between 0.5 and 0.7 /rHz. Only low degree modes 1=0,1 
and 2 have been taken into account in the computations. 


Composition 

Mass 

Age 

Teff 

m v 

Av 


(Mq) 

(Myr) 

(K) 


(pHz) 

Cl 

0.769 

625 

4707 

10.19 

206.2 

Cl 

0.8 

625 

4835 

9.92 

199.1 

Cl 

0.95 

625 

5454 

8.84 

170.2 

C2 

0.95 

625 

5375 

8.93 

171.5 

C3 

0.95 

625 

5297 

9.02 

172.7 

C4 

0.95 

625 

5612 

8.70 

167.0 

C5 

0.95 

625 

5232 

9.10 

173.1 

Cl 

0.95 

500 

5451 

8.84 

170.8 

Cl 

0.95 

700 

5456 

8.83 

169.8 

Cl 

1.00 

625 

5637 

8.54 

160.8 


for frequencies between 3.5 and 4.5 mHz. The lower limit is 
approximately the frequency of the mode of radial order n=15 
and degree f: = 0 in a 0.77 Mq model. For higher mass stars this 
frequency corresponds to higher orders because of the decrease 
of the acoustic radius with mass. It therefore always lies above 
n= 10 which keeps the explored domain far from possible radia- 
tively induced oscillations (see Sect. |2J. The other reason why 
we utilize frequencies in this range is that, as can be seen in Fig. 
0 the large separation reaches the asymptotic regime above ~ 
3 mHz; thus its direct interpretation in terms of acoustic radius 
becomes all the more reliable. Finally we set an upper limit of 
4.5 mHz on the computation of the mean large separation in or¬ 
der to reduce the scatter in this mean. We preferentially select 
the models at 625 Myr which has recently be estimated to be 
the Hyades age (Perryman et al. 1998). However we also ex¬ 
tend the analysis to 500 and 700 Myr to check for the impact 
of evolution slightly outside the claimed age uncertainties for 
the cluster ~ 50 Myr. Calculations of the large separation of 
lower mass stars (0.95 Mq) show its value to be extremely ac¬ 
curate. In the frequency range considered the large separations 
are extremely flat; the scatter is ~ 0.5/rHz (Fig. 0. Moreover 
we find that the large separation experiences very small varia¬ 
tions with age. The very slow evolution of stars having masses 
below a solar mass explains this stability. The 0.9 Mq star main 
sequence lifetime exceeds 14 Gyr (Schaller et al. 1992). We 
note that the large separation is also very stable with compo¬ 
sition (for any conceivable variation of the metal or helium 
fraction). We report a Av variation of ~ 5.7/rHz ie 3.4% of 
its total value for AT as high as 0.04. The different values for 
the helium fraction determined by Pinsonneault et al. (1998) 
(Y=0.283) and Perryman et al. (1998) (Y=0.26) lead to a Av 
variation of 3.2/rHz. The variation with metallicity is similarly 
small. Going from composition Cl to C5 there is roughly 0.2 
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Fig. 5. Large separation as function of frequency and for differ¬ 
ent masses; crosses : 0.95 Mq, diamonds : 1. Mq and triangles 
: 1.15Mq The vertical dashed lines shows the 3500-4500 //Hz 
range over which we consider the mean large separation. 

dex difference in metal content (iron group metals represent 
less than 8% of the global metal mass fraction). Then we note 
however that Av only changes by 2.9//Hz. It is well known that 
metals mostly affect the stellar plasma opacity. The opacity in 
turn changes the stellar stratification in radiative zones but has a 
much weaker effect in outer convection zones where stratifica¬ 
tion is mostly adiabatic. Yet this is precisely the part of the star 
that mainly contributes to a large separation value. As illus¬ 
trated by the aforementioned asymptotic expression the large 
separation is determined where the sound speed is at its lowest, 
i.e. in the coolest regions. In the 0.95 Mq star of Cl composi¬ 
tion the convection zone represents 0.29/?* but conversely its 
acoustic radius is 0.59 % of the total acoustic radius. For the 
other compositions we check that these values are different by 
less than 3 %. 

The age- and composition-related stability of Av make it 
easier for the effect of mass to show up. A change of 0.05 Mq 
gives rise to a « 1 ()// Hz discrepancy. This variation clearly is 
much larger than the variations associated with mass or compo¬ 
sition, and also larger than the resolution level of future space- 
borne experiments such as EDDINGTON that should deter¬ 
mine frequencies to the level of 0.3//Hz (Favata et al. 2000). 
Hence simple measurements of the large separation variation 
might be used for the mass determination of lower main se¬ 
quence single Hyades stars. The accuracy in mass determina¬ 
tion will reach a few percent of solar mass. As oscillations are 
expected for any solar-like star the sample of these objects hav¬ 
ing accurately determined masses will be greatly increased. 

The predicted mass accuracy is comparable to or above th 
usual precision achieved when one relies on binary systems. It 
improves towards lower mass objects where the age determi¬ 
nation uncertainties are negligible. However, as the mass de¬ 
creases some other difficulties arise because of higher com¬ 
plexity in the input equation of state or low-temperature opaci¬ 
ties (Chabrier & Baraffe 1997). Besides this many aspects of 
the physics of convection in the outer regions of cool stars 


(T e fF = 4000 - 5000 K) remain to be investigated (Ludwig et 
al. 1999; 2002). A straightforward interpretation of the large 
separation in terms of mass might seem optimistic and it is 
perharps safer to use the large separation as an indicator of 
the mass hierarchy. By chance the vB22 system provides the 
best known mass in the Hyades : its secondary component 
(0.769 ± 0.005 Mq; Peterson & Solensky 1988) lies in the re¬ 
gion we have studied (Tabled- 

4. The extension of outer convection in the Hyades 

The precise extension of convective mixing in stars is crucial in 
several respects. First it has a direct impact on stellar evolution. 
Namely the occurrence of lower overshooting would signifi¬ 
cantly affect the evolution of low mass objects after they leave 
the main sequence and when they ascend the red giant branch 
(Alongi et al. 1991). The history of the surface abundances of 
the light elements closely depends on the extent of the con¬ 
vection zone. This in turn strongly constrains any inference on 
the dynamical phenomena in envelopes (Piau & Turck-Chieze 
2002; Piau et al. 2003). Finally convection and rotation remain 
the key ingredients for improving our knowledge of the angu¬ 
lar momentum evolution or of the surface magnetic activity. 
In this instance observational studies of young stars are now 
capable of providing new clues on the precise nature of the dy¬ 
namo process responsible for surface magnetic manifestations 
(Feigelson et al. 2003; Donati et al. 2003). 

Thanks to the increase of computational power, direct 
three-dimensional modelling increasingly approches realistic 
physical conditions of stellar convection. Many of these studies 
have focused on solar-type outer convection zones (Brummell 
et al. 2002 and references therein) and now also begin to ad¬ 
dress convective cores (Browning et al. 2004). However, the 
corresponding convective flows are extremely swift and turbu¬ 
lent. Consequently it is not possible to model convection in its 
full complexity and a fortiori to model stellar evolution in the 
framework of such direct approaches. For almost half a cen¬ 
tury we have relied on parametrized theories for convection 
such as the mixing length theory (Bohm-Vitense 1958). As is 
well known the MLT possesses one main free parameter q/mlt 
that is empirically adjusted by means of measurements of the 
actual solar radius. Some studies of close binary systems sug¬ 
gest that no change of q/mlt is required from the Sun to other 
low mass stars (Fernandes et al. 1998) while others suggest 
some slight changes « 0.2 with mass, e.g. in the a Cen system 
(Eggenberger et al. 2004). There are no basic physical reasons 
why this parameter should be independent of age, mass or com¬ 
position. 

Our goal here is to check how seismic information can 
be used to acquire better knowledge of convection. On the 
one hand the convective efficiency in the superadiabatic lay¬ 
ers directly affects the extent of the convection zone. On the 
other hand seismology has provided direct measurements of 
the depth of the convection zone in the solar case (Christensen- 
Dalsgaard et al. 1991). It is therefore possible to relate seis¬ 
mic measurements to convective efficacy. Thanks to recently 
developed techniques one can address the issue of the depth 
of the outer convection zone in other stars than our Sun (e.g. 






L. Piau, J. Ballot and S. Turck-Chieze: Seismic constraints on open clusters 


9 


Monteiro et al. 2000; Roxburgh & Vorontsov 2001; Ballot et 
al. 2004, hereafter BTG04): the presence of a discontinuity in 
stellar structure is detectable through (low-degree) mode phase 
shifts where it induces oscillatory behaviour. For instance the 
second difference (cL= v^-i - 2 v^„ + n+ \) exhibits mod¬ 

ulations whose periodicities are directly related to the acous¬ 
tic depth of the discontinuities (see Fig. We quantitatively 
evaluate the variations of this parameter resulting from the con¬ 
vection / radiation discontinuity for various Hyades models in 
the 1 to 1.3 Mq range. This range lies above the masses ad¬ 
dressed in binary systems by Fernandes et al. (1998) and is 
interesting to address in the Hyades in several respects. As 
open clusters comprise stars with marginally different compo¬ 
sitions or ages, they offer the possibility to specifically check 
the impact of mass on the convective efficiency in the supera- 
diabatic layers. Then the oscillations will be more easily de¬ 
tected above ~ 1 Mq as the corresponding stars exihibit visual 
magnitudes below ~ 8 / 74 ' and show larger fluctuations than the 
Sun and their lower mass counterparts (cf Sect. Q. At the same 
time a Hyades member of ~ 1 Mq should have lost a signifi¬ 
cant fraction of its initial angular momentum and its rotation 
rate should be small all the way down to its center (Piau et al. 
2003). This is an important point as the present analysis does 
not take into account rotationally induced effects on models 
and oscillation frequencies. In this instance Hyades observa¬ 
tions suggest a rapid decrease of the equatorial velocity as a 
function of 77ir between 7000 and 6000 K. As suggested by 
Gaige (1993) whereas below 77fr = 6400 K there are no stars 
with an equatorial velocity above 30 km s _1 the stars above 
T e g = 7000 K exhibit a distribution with a large scatter from 0 
up to =3 200 km s -1 . For masses of 1, 1.15 and 1.3 Mq we com¬ 
pute r eff = 5640,6110,6500 K respectively at the age of the 
Hyades and the mass interval should therefore correspond to 
slow or moderate rotators. Finally the mass interval considered 
covers the red wing of the so called lithium gap (Boesgaard & 
Tripicco 1986): stars with the same age as the Hyades or older 
having effective temperatures between 6300 and 6850 K are 
lithium- and also beryllium-deficient in comparison to hotter 
or cooler stars. This is most probably related to a rotationally 
induced (magneto)hydrodynamic mixing process occurring at 
the top of the radiation region. For almost 20 years the lithium 
gap has been an active research field both theoretically (Talon 
& Charbonnel 1998) and observationally (Boesgaard & King 
2002 ). 

To determine the extent of the convection zone we make 
use of the second difference d^v^n, exclusively evaluated from 
modes of degree 7 = 0,1 and 2. For the order n we have con¬ 
sidered the range between 13 and 26 which corresponds to the 
most excited modes in the solar case. Two periods are clearly 
seen in 62 Vc, n vs vq„ (see Fig.|6j: the smallest oscillation is asso¬ 
ciated to the base of the convection zone (hereafter BCZ). The 
largest corresponds to the discontinuity induced by the second 
ionization region of helium which has been discussed in Sect. 
1 m To reduce its impact we subtract a smoothed interpolation 
of the raw data from 62 Vc, n - This procedure is more precisely 
described in BTG04. We then perform a spectral analysis to 
evaluate the acoustic depth of the convection zone tbcz (Fig. 
[6). The transform gives tbcz which is an estimate of tbcz bi¬ 
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Fig. 6. Upper panel : the second difference as function of fre¬ 
quency for two Hyades models (625 Myr) of similar 77 jr ~ 
6550 K but with different ctmlt and masses. <> :(- (), 
A:7=l,m:7’ = 2. Dashed line : 1.35 Mq and £Zmlt = 1.566. 
Continuous line : 1.3 Mq and czmlt = 1.966. In order to make 
the figure clearer the upper and lower tracks have been shifted 
by 1.5 and -1.5 pHz respectively along the ordinate axis with 
respect to the original values. Lower panel: amplitudes for the 
Fourier transforms of second differences presented in the upper 
panel with similar conventions on linestyle. A smooth interpo¬ 
lation of second differences has been performed and subtracted 
to raw data before the Fourier transform in order to attenuate 
the He II ionization effect (this effect is still visible through the 
peaks around 800 s). The vertical dashed line shows the middle 
of the continuous line peak at half maximum. 


ased by surface effects (Christensen-Dalsgaard et al. 1995). In 
the solar case there is for instance a 200 s difference between 
tbcz and tbcz- Similarly half the inverse of the large separa¬ 
tion, is shifted by this amount with respect to the acoustic 
radius tr (cf equation|6}. The difference -f B cz is therefore 
unbiased and corresponds to the acoustic radius of the convec¬ 
tive base At bcz (BTG04). 


Atbcz - tr - tbcz — tt. - tbcz 

2Av 


(7) 
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Fig. 7. Acoustic radii of the bases of the outer convection zones 
scaled by Av/Avq and as a function of the effective tempera¬ 
ture at the age of the Hyades. Three qtmlt have been consid¬ 
ered.The symbols show results extracted from 62 V. Triangles : 
a M LT-l-566, masses 0.95,1,1.05,1.15,1.3 and 1.35 Mq; Stars 
: qtmlt- 1-766, masses 0.95, 1, 1.15 and 1.3 Mq. Diamonds : 
<zmlt= 1-966, masses 0.95 , 1, 1.15 and 1.3 Mq. The dashed 
lines indicate the true BCZ positions in the different models 
The horizontal error bar cr(T e g) is 2 x 55 K. The vertical error 
bars are discussed in the text. 

We use A t bcz as it is not (or less) dependent than f BCZ on 
surface effects where the physical description is less under con¬ 
trol (turbulent pressure, superadiatic effects). Figure |7] shows 
the acoustic radii of convective bases A t B cz scaled by Av/Avq 
as a function of effective temperatures and for the three o'mlt 
values 1.566, 1.766 (our solar calibrated value) and 1.966. The 
acoustic radii scales with the mean density and inversely with 
frequencies. As in Fig.0 the scaling factor Av/Avq therefore 
makes it possible to isolate the effect of the o'mlt parame¬ 
ter. It moreover reinforces our ability to determine changes in 
A t B cz with effective temperature because it extends the varia¬ 
tion range of A t B cz- The higher o'mlt higher the convective ef¬ 
ficiency and deeper the convection zone. As observed on Fig.0 
this in turn translates into smaller acoustic radii for the internal 
radiative regions. On Fig. 0 we plot a typical cr(T e ff) error bar 
of 55 K for the effective temperature as estimated by Thorburn 
et al. (1993). The errors made in the determination of A t B cz 
are extensively discussed in BTG04. We summarize the results 
here. Firstly there is a intrinsic error due to the method itself. 
This includes the effects of the limited number of points and the 
unknown residual bias of the estimator. For the modelled solar¬ 
like stars, this error is smaller than 1.5%, and generally smaller 
than 1 % in units of total stellar acoustic radius. This difference 
can be observed in Fig. □ by comparing the symbol positions 
with respect to the dashed lines. Secondly the effects of the 
uncertainties on the frequencies can be estimated with Monte 
Carlo simulations. For expected errors of 0.1-0.15 pHz on the 
frequencies, the error induced in the determination of A t B cz is 
around 2-2.5 %. Thus we have considered a pessimistic global 
error of 3 % in units of total stellar acoustic radius. 

Using a 2D hydrodynamic code including radiative transfer 
Ludwig et al. (1999) evaluate the stratification of the outer con¬ 


vection zone for various surface conditions in T e ff and log g. 
They cover stellar types from K to F (4300K < T e g < 7100 K 
and 2.54 < log g < 4.74). The o'mlt should vary from ~ 1.3 to 
» 1.75 to match the entropy jumps they find. 

The 0.4 range we consider here for qtmlt is therefore repre¬ 
sentative of the variations with surface conditions expected in 
convection efficiency. Interestingly Lebreton et al. (2001) find 
that Hyades isochrones having o'mlt = 1-6 - 1.8 (respectively 
o'mlt ^ 1-4) better fit the My - B - V observed magnitude 
diagram in the 1.3 - 1.7Mq range (respectively below IMq). 
We note here that owing to the errors any variation of convec¬ 
tive efficiency with stellar type corresponding to an amplitude 
Ao -mlt ~ 0.2 should be clearly observable in a diagram similar 
to Fig0 

To push the investigation further we have computed stellar 
models with variable o'mlt parameters. Following the Ludwig 
et al. (1999) scaling-law, o'mlt varies with the surface condi¬ 
tions (T e ff and log g). With regard to constant (solar-calibrated) 
o'mlt = 1-76 models the difference in Atbcz for the variable 
o'mlt models is quite small (5 20 s) except for the 1.3 Mq star 
where it reaches 95 s. This stems from Ludwig’s result that 
around solar conditions the o'mlt values exhibit a plateau in 
the T eS - log g plane. The slightly larger Atbcz (~ 100 s) of 
the ‘variable 1 o'mlt L3 Mq star reflects the decrease in convec¬ 
tive efficiency towards higher T e g predicted by Ludwig et al. 
(1999). Let us stress in this example the ability of seismology 
to give constraints for the hydrodynamic convection simula¬ 
tions (Robinson et al. 2004). 

Table^summarizes some quantities characterising the con¬ 
vection together with the corresponding stellar surface condi¬ 
tions (T eff and log g). We recall here that although o'mlt surely 
is the most important parameter within the MLT framework it 
is by no means the only one. There are indeed several formu¬ 
lations of the MLT that give different stratification for identical 
<*mlt values as illustrated by Ludwig et al. (2002). To specify 
the MLT formulation we use, we provide four other parameters 
(such as the surface to volume ratio for convective elements) in 
a short appendix. Like any other local approach of convection 
(eg Canuto & Mazzitelli 1991), the MLT is simplistic and can 
be criticized on many points. Indeed various o'mlt parameters 
are required to match hydrodynamic computations for velocity 
or thermal stratification in the upper convection region (Ludwig 
et al. 2002). For these reasons, we provide a more meaningful 
physical quantity than o'mlt : the specific entropy of the deep 
convection zone s cz . By entropy of the deep convection zone 
we mean the entropy level reached where convective motions 
become extremely close to adiabatic. s cz is directly related to 
the extent of the convection zone and to the convective effi¬ 
ciency in superadiabatic layers. It may therefore be assigned 
an intermediate role between hydrodynamic computations and 
seismology. 

We shall end this section with two caveats. Firstly and as 
noticed in the introduction no overshoot was considered here. 
Monteiro et al. (2000) show that such a phenomenon would in¬ 
crease the amplitude and the acoustic depth of the signal result¬ 
ing from the BCZ discontinuity. Thus any interpretation of the 
extent of the outer convection zone in terms of convective effi¬ 
ciency or o'mlt parameter depends on the overshoot. We never- 
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Table 4. Column 1 : mass. Columns 2 & 3 : surface effective temperature T e g and gravity log g. Column 4 : o-mlt, column 5 : 
acoustic radius of the base of the outer convection zone A tbcz as deduced from the Fourier transform of the second difference 
and typical error bars; column 6 : specific entropy in the deep convective region s cz . The models constructed with a variable ctmct 
are to be found in the last section of the table. The reported »mlt values are those corresponding to the surface conditions at the 
age of the Hyades. 


Mass 

Teff 

log g 

®MLT 

A t BC z 

s cz 

(Mq) 

(K) 



(s) 

(erg K 1 g-'/lO 9 ) 

0.95 

5386 

4.54 

1.566 

1155±122 

1.67 

1 . 

5561 

4.51 

1.566 

1304±120 

1.72 

1.05 

5725 

4.48 

1.566 

1426±152 

1.77 

1.15 

6025 

4.41 

1.566 

1892±157 

1.88 

1.3 

6426 

4.31 

1.566 

2934±198 

2.19 

1.35 

6545 

4.28 

1.566 

3394±194 

2.38 

0.95 

5454 

4.56 

1.766 

1123± 121 

1.66 

1 . 

5637 

4.53 

1.766 

1232±118 

1.70 

1.15 

6111 

4.44 

1.766 

1777±165 

1.85 

1.3 

6504 

4.34 

1.766 

2697±206 

2.12 

0.95 

5511 

4.57 

1.966 

1085±120 

1.65 

1 . 

5702 

4.55 

1.966 

1185±119 

1.69 

1.15 

6186 

4.46 

1.966 

1609±130 

1.83 

1.3 

6575 

4.35 

1.966 

2419±216 

2.39 

0.95 

5458 

4.56 

1.781 

1120±121 

1.66 

1 . 

5640 

4.53 

1.774 

1209±130 

1.70 

1.05 

5807 

4.50 

1.774 

1387±142 

1.74 

1.15 

6101 

4.43 

1.741 

1769±164 

1.86 

1.3 

6465 

4.32 

1.662 

2792±193 

2.15 


theless note that 1) the actual helioseismic inferences on over¬ 
shooting suggest that this phenomenon is moderate, Monteiro 
et al. (1994) find l ov = 0.07 H v ; 2) hydrodynamic simulations 
of the solar convection zone similarly suggest rather moderate 
overshooting : below 0.1 \H V (Brummell et al. 2002). The pre¬ 
dicted movements moreover do not mix the stable region with 
enough efficiency to change its stratification significantly 7 * i.e. 
leave no seismic signature. Secondly the diffusion process has 
been taken into account following Michaud & Proffitt (1993). 
Although we do not report it here in details we mention that the 
amplitude of the BCZ-induced oscillatory behaviour in So vy.„ 
increases both with age and mass of our Hyades models. This 
is opposite to the comparable effects on frequencies reported 
in Monteiro et al. (2000). Microscopic diffusion is responsi¬ 
ble for this difference. Monteiro et al. mostly addressed ZAMS 
stars where this process has no impact and was therefore not 
included in their computations. The description of diffusion we 
use is valid for G dwarfs and later types but becomes unrealistic 
for earlier types. Above T e ff ~ 6500 K the radiative diffusivity 
which we do not consider plays a non-negligible role (Morel & 
Thevenin 2002). Besides this the stellar rotation increases and 
the microscopic diffusion interacts more strongly with one or 


7 they represent therefore strictly speaking overshoot and not pene¬ 

trative convection. 


several macroscopic mixing processes as it is suggested by the 
occurrence of the lithium dip. 

5. Summary and conclusion 

We computed classical and seismic stellar parameters using 
the stellar evolution code CESAM and the adiabatic pulsa¬ 
tion package of Christensen-Dalsgaard. Classical and seismic 
variables were employed to constraint the stellar models. Only 
{ = 0,1,2 low-degree modes were considered. As a prelimi¬ 
nary analysis we estimated apparent visual magnitudes, lumi¬ 
nosity fluctuations and a few other interesting seismic variables 
in three nearby open clusters. We subsequently focused on stars 
within the closest and best-known of these : the Hyades. We 
predicted (Sect. □ that the Hyades stars with masses around or 
above solar should exhibit oscillations detectable by the near¬ 
future seismic space-borne experiments. Future ground-based 
observations with HARPS or ESPaDOnS, in probing surface 
velocities could similarly begin to answer some of the ques¬ 
tions mentioned here. Exploiting the predicted oscillation fre¬ 
quencies we had three goals: to set constraints on mass, compo¬ 
sition and structure. Let us briefly sum up our results on these 
points and conclude. 

The mean value of the large separation Av in the 3.5 to 4.5 
mHz frequency range is a strong function of the mass: we cal- 
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culated (Sect. that a « 0.05Mq mass change induces a 
~ 1 0 /j Hz variation in Av. This is much higher than the Av vari¬ 
ation expected from typical uncertainties in age or composition 
(Table Asteroseismology will therefore allow the determi¬ 
nation of absolute and relative masses in nearby open cluster 
and with competitive accuracy when compared to the best as¬ 
trometric and spectrometric mass evaluations to date. However, 
at variance with these measurements seismology requires only 
short observation runs and is not limited to a very small sample 
of binaries. 

Asteroseismology can shed new light on the composition 
of open clusters(Sect. ED . As main opacity contributors, met¬ 
als strongly influence the effective temperature T eff and have 
a direct impact on the internal structure which affects the large 
separation Av. Using Hyades models of different masses we 
have shown that the location of the T eff - Av relation in the 
corresponding diagram provides a measurement of the global 
metal content. In particular, we are able to distinguish models 
buildt with solar-intermetallic ratios from models buildt with 
non-solar-intermetallic ratios. This is an important issue as we 
recall that 1) iron only represents « 2 % of the total metal con¬ 
tent in number 2) there is a clear scatter in [O/Fe] vs [Fe/H] in 
the Galactic disk (Edvardsson et al 1993). As is well known, 
the amplitude A He of the second difference Sov main modula¬ 
tion stems from the second ionization region of helium. Its use 
to infer helium fraction in the envelope however requires other 
stellar parameters such as mass and radius to be known (Basu 
et al. 2004). The FINSEN 342 star belongs to a Hyades binary 
system. Its mass and age are accurately known. Given these 
constraints, we showed that for this star a variation of the ini¬ 
tial helium mass fraction from 0.26 to 0.28 would increase An e 
by 0.2/rHz. 

The main feature of structure of solar-like stars is the pres¬ 
ence of a radiative core surrounded by a convective region. 
Using 62 V modulations with frequency we have evaluated the 
acoustic radii of the bases of the convection zones (BCZ) for 1 
to 1.3 Mq Hyades models and to an accuracy better than 3 % of 
the total acoustic radius (Sect. 0. As, for a given effective tem¬ 
perature, the acoustic radius of the BCZ depends on convective 
efficiency, asteroseismology provides a possible tool for con¬ 
straining a M LT (Fig. m As is well known, the mixing-length 
theory does not give a physically consistent picture of convec¬ 
tion. Nevertheless, it is interesting to constrain otmlt in stars as 
the ongoing hydrodynamic computations on convection make 
predictions for this parameter. Seismology therefore offers a 
means to check our calculations (and understanding) of con¬ 
vection in the superadiabatic layers. The considered mass range 
is interesting in several respects : it starts where the knowledge 
of our Sun can presumably be transposed and then extends in a 
domain where variations in convection efficiency are suspected 
through HR diagram analysis (Lebreton et al. 2001) or binary 
systems calibrations (Eggenberger et al. 2004). Moreover the 
interval covers the red wing of the lithium dip where additional 
mixing is at work in the radiation zone (Talon & Charbonnel 
1998). Setting constraints on convection depth in this range 
will not only shed new light on the convection efficacy in re¬ 
lated stars but will also help constraining any deeper dynamical 
mixing phenomenon. 


Because they comprise stars of similar age and composi¬ 
tion, open cluster are unique laboratories for stellar physics. In 
this respect the Hyades have proved to be of primary impor¬ 
tance. This importance is transposed to the asteroseismology 
field. In the Hyades it is possible to take advantage of collective 
effects and also to select some especially interesting binaries. 
In this work we presented various seismic tests that could be 
performed assuming today’s accuracy on the measurement of 
oscillation frequencies. These tests clearly open up new oppor¬ 
tunities in the field of stellar physics. Applied to the Hyades, 
asteroseismology will allow precise mass and composition de¬ 
terminations which are crucial for valuable astrophysical tests. 
It will furthermore enable measurements of the depths of outer 
convection zones which in turn will provide indications of con¬ 
vective efficiency in stars. 
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Appendix 


The MLT version we use relies on the following equations (Cox & 
Giuli 1968) where fa, fa, fa and fa are dimensionless free parameters 
each having a specific physical meaning : 

2 _ s A 2 git(v-v f ) 

V Hf 

v is the convection velocity, fa describes the dissipation of the 
work of buoyancy forces due to viscosity and we set it to . 


p _ r pcpi/TA(V-V e ) 

t ~ J 1 Hp 

F is the convective energy flux. / 2 is related to the proportion 
between upward and downward matter movements and is set to -1. 


r = v7 
L is 


_ P C P 


ad 

the 


r t -y u ‘ e /1 , i4\ 

/3<rr- ,u + r 2 

convection efficiency. 


fa is related to the 
surface I A. volume ratio of convective elements and is set to 
24. Finally fa is related to radiative transfer inside convective bubbles 
and is set to 3. Most of the symbols in the preceding equations keep 
their usual meanings (eg A = qtmltHp is the mixing length) and 
were defined in Sect. 0 In addition : T is the temperature; a is the 
Stefan-Boltzmann constant; V, and V„d are the mean thermal 
gradient, the convective element thermal gradient and the adiabatic 
thermal gradient respectively; r e is the Rosseland optical thickness 
of a convective element. The four dimensionless parameters are 
similar to the Bohm-Vitense (1958) formalism with the exception 
of fa which she assessed to be zero. As remarked by Henyey et al. 
(1965), however, the value of this parameter only has a marginal 
impact on the general description of convection through MLT. We 
have moreover directly checked that our approach and results indeed 
are similar to other computations done within the Bohm-Vitense 
formalism. 
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